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Abstract 

We investigate the sub-leading contributions to the free energy of Bethe Ansatz solv- 
able (continuum) models with different boundary conditions. We show that the Ther- 
modynamic Bethe Ansatz approach is capable of providing the 0(1) pieces if both the 
density of states in rapidity space and the quadratic fluctuations around the saddle 
point solution to the TBA are properly taken into account. In relativistic bound- 
ary QFT the 0{1) contributions are directly related to the exact (/-function. In this 
paper we provide an all-orders proof of the previous results of P. Dorey et al. on 
the (^-function in both massive and massless models. In addition, we derive a new 
result for the (^-function which applies to massless theories with arbitrary diagonal 
scattering in the bulk. 
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1 Introduction 

The study of the thermodynamics of one dimensional integrable models with factorized scat- 
tering dates back to the seminal work of Yang and Yang [1, 2]. Their method, known today as 
the Thermodynamic Bethe Ansatz (TBA) is quite general and it was worked out for a large 
number of models relevant to condensed matter physics [3]. In its simplest formulation the TBA 
is written down for periodic boundary conditions and it provides the free energy density, ie. the 
0{L) part of the free energy. 

In this paper we study the sub-leading pieces of the free energy for different boundary condi- 
tions in continuum models. We restrict ourselves to theories with diagonal scattering; however it 
is expected, that the proposed methods will work even in the non-diagonal case ^. In this paper 
we focus mainly on relativistic theories. 

In integrable relativistic field theory the TBA was introduced by Al. Zamolodchikov in [4]; 
soon thereafter it became one of the central tools to study finite size effects. In relativistic models 
Euclidean invariance implies that the free energy density at finite temperature is directly related 
to the exact ground state energy in finite volume. Studying this quantity it is possible to recover 
the behaviour around the fixed points of the renormalization group flow, which are usually given 
by a conformal field theory. Thus the TBA provides a link between the scattering theory (IR) and 
perturbed CFT (UV) description of the same model: it predicts the central charge, the scaling 
dimensions of the perturbing operator, and various other quantities [4, 5, 6]. 

^see the related comments in the Conclusions (Sec. 6) 
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The techniques of integrabihty can be apphed to problems with non-trivial integrable bound- 
ary conditions; in relativistic scattering theory the foundations were laid down in [7]. One object 
of particular interest is the exact g-function, which is the off-critical generalization of the non- 
integer ground state degeneracy of critical boundary conditions introduced by Affleck and Ludwig 
in the context of the Kondo model [8]. The (^-function describes the 0{1) contribution of a sin- 
gle boundary to the free energy and it can be used to study renormalization group flows in 
the space of boundary field theories [9, 10]. In [11] it was shown that the (7-function satisfies 
a gradient formula, from which it follows that in unitary theories the boundary entropy mono- 
tonically decreases under the RG fiow. This is the ^(-theorem, which can be regarded as the 
boundary-counterpart of the celebrated c-theorem by A. B. Zamolodchikov [12]. 

It is an old and very natural idea to determine the g-function in the framework of the 
Thermodynamic Bethe Ansatz. The first result appeared in [13], where the authors proposed a 
simple formula based on the boundary-dependence of the Bethe equations. Later it was found 
in [10] that although the results of [13] correctly describe the boundary-dependence of the g- 
function, a boundary-independent term has to be added in order to match the predictions of CFT. 
The missing piece was derived in [14] using a cluster expansion for the free energy; the exact result 
was expressed in terms of the solution of the TBA with periodic boundary conditions. While this 
exact 5f-function successfully passed a number of non-trivial tests [14, 15] and recently it was 
generalized to describe a massless flow in [16], the interpretation of the boundary- independent 
terms remained unclear. 

A remarkable attempt to obtain the non-extensive pieces to the free energy of Bethe Ansatz 
systems was performed in [17], where it was shown that the quadratic fluctuations around the 
saddle point solution yield a well-deflned 0(1) piece. However, the calculation of [17] seemed 
to contradict all previous results: it did not reproduce the boundary-independent term of [14], 
moreover it predicted an 0(1) piece even in the periodic case where no such term is expected. 

In this paper we revisit the calculations of [17] and argue that the only flaw of [17] is that 
it did not take into account the non-trivial density of states in the conflguration space of Bethe 
Ansatz systems. In other words, the functional integral for the partition function was built on 
an incorrect integration measure. We propose a new normalization based on the thermodynamic 
behaviour of the density of states and we obtain the correct results in all previously considered 
cases. 

The paper is organized as follows. In the next subsection we provide the necessary deflnitions 
for the g-function of relativistic boundary field theory. Sec. 2 serves as a warm-up: we consider 
general Bethe Ansatz systems and the behaviour of the density of states in the thermodynamic 
limit. In. 2.2 we revisit the calculations of [13] and show that it is possible to obtain the boundary- 
independent part of the (7-function by a simple heuristic argument. Motivated by these findings 
in Sec. 3 we present a general framework to evaluate all 0(1) pieces to the free energy. These 
formal results are then evaluated explicitly in massive and massless relativistic models in Sections 
4 and 5, respectively. Finally Sec. 6 includes our conclusions. 

1.1 The exact (51-function — definitions 

The exact g(-function can be defined as follows [13, 14]. We restrict ourselves to the sim- 
plest case with only one massive particle in the spectrum; the generalization to other models is 
straightforward. 

Let us consider a finite cylinder with height L and circumference R (fig. 1). The integrable 
boundaries a and b are placed on the two ends of the cylinder. The partition function can be 
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Figure 1: Pictorial representation of the two channels for the evaluation of the cylinder partition 
function. The (imaginary) time evolution is generated in the vertical direction by the correspond- 
ing Hamiltonians. 



evaluated in two different channels. Viewing R as the direction of time (fig. 1(a)) one obtains 

ZabiR, L) = Tr e-^-(^)« = J] e'^'^^)^, (1.1) 

where Hab{L) is the Hamiltonian of the open system of size L with integrable boundary conditions 
described by the reflection factors Ra{0) and Rb{0). The summation runs over a complete set of 
states and E'^{L) are the eigenvalues of the Hamiltonian Hab{L). 

On the other hand, one can use L as the time variable (fig. 1(b)). In this picture the bound- 
aries play the role of initial and final states of the time-evolution operator H(R), which is the 
Hamiltonian of the system of size R with periodic boundary conditions. In this channel the 
partition function is evaluated as 

Zab{R,L) = {Ba\e~^(^^^\Bb) = {Gt{R)yGt{R)e-^^^^^^ (1.2) 

Here \Ba) and \Bb) are the boundary states corresponding to the boundary conditions [7] and 
the summation runs over a complete set of states of the periodic be. system. The amplitudes 
G*j(i?) are defined as the normalized overlaps 

Gt{R) = Wl, J = -^b (1.3) 



Vim' 



Equations (1.1) and (1.2) have to be contrasted with the definition of the partition function with 
periodic boundary conditions in both direction: 

Z{R, L) = Tr e-^(«)^ = e'^''^^)^ = Z{L, R) (1.4) 

If mL ^ 1 then (1.2) and (1.4) are dominated by the ground state with energy 

Eo{R) = eR + 0{e-'^'^), 
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where e is the bulk energy density, which can be determined by comparing the TBA results to 
conformal perturbation theory [4, 5, 6]. 

Similarly, if mR S> 1 then (1.1) is dominated by the ground state of the boundary system 
with energy 

<(L) = 6L + A + /, + 0(e— 

where fa and /f, are non-extensive boundary-contributions which can be obtained from the bound- 
ary TBA [13, 9]. 

Comparing (1.1) and (1.2) in the regime mL,mR ^ 1 one finds 

G°.(R) = e-/j^(l + 0(e-"^^)) , j = a,b 

The 5f-function is then traditionally defined as 

GO(i?) = e-^^%(i?), j = a,b 

It follows from dimensional arguments that the (/-function depends on r = mR only. It is useful 
to re-define the partition functions as 

Z(i?,L) = J]e-(^v.W-^)^ (1.5) 

With this prescription the vacuum energies have the asymptotics 

lim E^^{L) = lim Eq{L) = 0, 

L—^oo L-^oD 

and the excited state energies are calculated additively in the Bethe Ansatz picture. The g- 
function is then given by the limit 

log ga{r)gb{r)= lim (log Zab{R, L) - log Z{R, L)) (1.7) 

Our goal is the evaluation of the partition functions (1.5) and (1.6) in the framework of the 
Thermodynamic Bethe Ansatz, with a special emphasize on the overall normalization. The exact 
g'-function will be determined by the relation above. 

We would like to remind the reader that in relativistic scattering theories the Bethe Ansatz 
does not provide an exact description of the spectrum. In fact there are exponentially decaying 
residual finite size effects, which modify the multi-particle energies obtained from the Bethe 
Ansatz [18, 19]. In the conventional derivation of the TBA these contributions are neglected; 
this is certainly a good approximation in the dilute regime, where the average distance between 
the particles is much larger than their Compton-wavelengths. However, it can be argued that 
the TBA yields the correct results at any temperatures and any densities [20]. In this paper 
we develop arguments about the thermodynamic behaviour of certain Bethe Ansatz quantities. 
We do not use any special assumptions other than those already used in the derivation of the 
TBA. Therefore we may neglect exponential corrections; this assumption will be justified by the 
results, which are in agreement with the CFT results in the UV (high temperature) limit. 
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2 Thermodynamic Bethe Ansatz and density of states 

Let us consider Bethe Ansatz systems with one particle type and no internal degrees of 
freedom. We use the rapidity variable to parametrize the states. The scattering is assumed to 
be elastic and factorizing; the two-particle scattering is described by the pure phase 

S{ei -9j) = e'^^^^-^^^ 

Energy and momentum are given by the functions e{9) and p{0). In relativistic field theories 

e{9) = ?ncosh0 p{9) = msinh^ 
whereas in the non-relativistic case we have 

e{e) = p{6) = mO 

Alternatively, one can introduce a chemical potential according to 

e{e) e{e) - n 

In the relativistic case we always set /i = 0. 

Consider the Bethe- Yang quantization of an A^-particle state in finite volume L with periodic 
boundary conditions: 

-Q g^g^ -0,) = l, j = l... N, 

where pj = p{9j). In the logarithmic form: 

Qj =PjL + ^^{ej -Ok) = 2TrIj j = l...N, Ij £ Z (2.1) 

We assume that the quantum numbers {/i, . . . ,Ii\f} completely characterise the state and that 
the solutions consist of purely real rapidities, ie. we do not consider theories with string- like 
solutions. Moreover, we assume that the solutions of (2.1) span the whole Hilbert-space. 
The multi-particle energies can be calculated additively: 

N 

E[{h,...,lN}]=^e{ej) 
i=i 

The solutions of (2.1) are placed evenly in the space of the quantum numbers. In rapidity 
space one can define the density of states as 

Pn{Oi, . . . ,6n) = det J, Jik 



Alternatively, this can be interpreted as the norm of the Bethe Ansatz state [21]. 

In the presence of integrable boundaries the quantization conditions can be written as 

e^2f.ii?„(^.)i?,(0^.) J] 5(% - ek)S{e, + 0fe) = 1, j = 1 . . . iV, > (2.2) 



or 



Q,- = 2p.jL + ^a{Oj) + MO,) + - Ok) + H0j + Ok) = 27rl,, j = l...N, 6, > 

k+3 



6 



where iia,fe(^) = e^"^"-'^^^^ are the elastic reflection factors of the boundaries. The density of states 
is defined as 

PAr(6'i, . . . ,0Ar) = det jf, Jik = -^ 

We will be interested in the thermodynamic behaviour of pjy a-^d pj^j. These quantities will 
play an important role in the normalization of the partition function in 2.2 and in the general 
treatment of section 3. 

If the particle number is kept fixed and the volume is sent to infinity, then the densities 
behave as . In particular 

TV 



PN,PN ~ + X J](27ra(^j)L 



where (j{0) = p'(0)/27r. Note that although the leading piece is the same for p]\f and p^v, the 
sub-leading terms are in general different. An important property is that the coefficient of these 
0{L~^) term grows linearly with N . Therefore it is expected that if both L and go to infinity 
with their ratios fixed, then the densities will have a non-trivial thermodynamic limit, which is 
expected to be different for /jtv and p^. In this limit all higher order terms of 0{L~'^) contribute. 

In the following we consider the particular case when the state in question is a finite- volume 
micro-canonical realization of an infinite volume thermal state. In other words, we send and L 
to infinity with their ratios fixed, and we choose the distribution of Bethe-Ansatz roots according 
to some Thermodynamic Bethe Ansatz equation. 

As a starting point we briefly review the basic notions of the TBA. In the thermodynamic 
limit we introduce the total density of states^ p{9) and also the density of occupied states p^°\9). 
The particle density is given by the integral 

n = ^= I de p^°\e) 



L 

In the periodic case the densities satisfy the constraint 

de' 



p^e) = -ap{e) + j ^ —^{e - e')p^°){e'\ (2.3) 

where (p = —i-^ log 5(0) and ap{e) = dp/de. On the other hand, in the boundary case we have 
e > and the densities satisfy the constraint 

1 /"OO I/D/ 

p{e) = -cjam + —{^{e - e') + ^{e + e'))p^°\e') (2.4) 

where now 

'"'^'^^^ " le (^^^^^ " 4 ^"'^^0 " ^""^^^^ ^^-^^ 

where Rab{e) = Ra{e)Rh{e)S{—2e) and the extra term containing 6{e) was introduced to cancel 
the unphysical solutions corresponding to = 0. 



^We use the term "density" for both pjv and p{6), p^°\B). However, these quantities have a very different 
meaning: p{8) and po{8) describe the distribution of roots within one configuration, whereas pjv describes the 
density of states in the (thermodynamic dimensional) configuration space. 
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It is known, that the the distribution of roots is given in both cases by the thermodynamic 
Bethe Ansatz equations 

/OO JQl 
^_^(0_^')log(i + e-(^')), (2.6) 



where we introduced the pseudo-energy as 

P{0) 



1 + e'^^'^ (2.7) 



2.1 The thermodynamic hmit of the densities and 

The calculations will be based on the techniques developed in the Algebraic Bethe Ansatz 
literature, see for example Appendix A in [22]. 

In the periodic case the matrix elements of J read 

N 

Jik = 5ik{Lap{ei) + ^(/^(^i - 0j)) - '^{Oi - 9j) 
i=i 

The matrix J' can be written as the product 

J = Ge, where 

Qij = Sij-fj, Gij = 5ij - f^^^ and 

N 

^,=Lap{e,) + Y,v{e,,) (2.8) 

i=l 

With this notation 

Pn = det G N det 6 at 
In the L — 7- OO limit we have from (2.3) 

7j 27r Lp{9j) 

The elements of Gn can be written asymptotically as 



Gii — 5, 



1 (^(%) 



I] — ^ij 



2-kL p{9j) 

Using (2.7) we conclude that the limit of det Gij is given by the Fredholm determinant 

det (i - P^) (2.9) 
where P~ is an integral operator acting as 

Therefore 

N 

PN =^ det (^i - X JJ 27rL/3(6'j) (2.11) 
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The Fredholm determinant (2.9) is well-defined if the operator P~ is trace class. This can be 
checked explicitly in non-relativistic situations [21], or in massive relativistic models at low 
temperatures, where 

Trf (p-yp-) ~ 0(e-2-/^). 



Here and in the rest of the paper we will assume that the Fredholm determinants we encounter 
are always well-defined. Even if this is not the case, they could be made regular by introducing 
a rapidity cut-off. This cut-off would not affect the final results, because the pieces of the 

free energy will be expressed as integral series which are well-defined for arbitrary temperatures. 

Now we turn to the calculation of p^r. In this case the matrix elements of the Jacobian read 

N 

i=i 

The matrix can be written as 

J = Ge, with 

i'3 — U 13 ' *J — *J ' 

7i 



where now 



N 



7,- = LaabiOj) + {^{e^ - 9,) + ^{Oi + 9,)) (2.12) 



i=l 

The quantities 7^- have the same thermodynamic limit as in the periodic case: 

7, ^ 2^Lp{9,) 

In the thermodynamic limit one finds 

N 



PN 



det {i - Q~) X JJ 2TTLp{9j), (2.13) 



where Q acts on functions defined on M"^ as 



(Q-(/))(x) = l^^i^i^ -y)- ^(^ + y)) (2-14) 

Notice that the leading parts of (2.11) and (2.13) have the same form, however the pre-factors 
are different. This will play an important role in the next subsection. 

2.2 A heuristic derivation of the function (massive case) 

Here we revisit the calculation of [13] to determine the exact ^r- function in massive relativistic 
boundary QFT. We obtain obtain the correct boundary independent part of the (7-function using 
a simple heuristic argument. 

In section 4 of [13] the authors start with the quantization conditions for N particles ^ 

e'^'^''''^'''^Ra{9,)R,{9,) n S{9, - eu)S{9, + 0^) = 1, i = 1 . . . iV, 0, > (2.15) 



^In [13] (and also in [14]) the roles of L and R are switched as opposed to our conventions 
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This can be written alternatively as a quantization condition for 2n particles: 

"'^"'^'''''^^^R{ei)Y[S{ej -Ok) = h j = l...2n (2.16) 



with the additional constraint 

The function R{0) in (2.16) is defined as 



(2.17) 



R{e) = Ra{e)Rk{e)s{-2e), 

In [13] the authors use (2.16) as a starting point to derive the thermodynamics. They show 
that (2.16) yields the usual periodic-boundary-conditions TBA and derive a single 0(1) correction 
to the free energy: 

1 1'^ 

\og{gagb) = ^J de @ab{0) log (1 + e-^(^)) (Incomplete) (2.18) 

where 

Qab{e) = -i^\ogR{e)-27:5{e) 

where the term 5{6) is introduced to cancel the unphysical states which do not respect the 
Pauli-principle. This result was shown to be incomplete in [10, 14]. 

We believe that the derivation in [13] is correct in every respect but one. There it was assumed 
that the change from system (2.15) to (2.16)-(2.17) is trivial. However, the previous subsection 
shows that the density of states of the two BA systems behaves differently in the thermodynamic 
limit. In fact, the periodic system with 2A^ particles has a density 



2N 

P2N ~ det (^i - P-) X J]^ 2^Lp{ej, 



whereas the two copies of the boundary systems with N particles have a density 



PNj ~ det^ - Q-j X JJ 2TiLp{9 



In other words, there is a finite ratio between the (thermodynamic dimensional) densities: 

P'^N det f i - P- 



2.19) 



This factor does not grow with the volume L, therefore it does not affect the thermodynamic 
limit of the distribution of roots. This means that passing from the system (2.15) to (2.16)-(2.17) 
the saddle point will not be changed and the distribution of roots will be described indeed by 
the periodic-boundary-conditions TBA. However, a finite factor for the density of states has to 
be kept in order to have a direct comparison between the partition functions. In other words, 
the ratio (2.19) has to be included in the overall normalization of the partition function of the 
boundary system. 
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It is useful to explicitly evaluate this finite factor. We start by employing the identity 




det ( i - K I = exp <^ - y -Tri^" } (2.20) 



which is valid for any Fredholm-determinant. In the present case we have 

^2 



n.A^ det^ 1-Q- r °° 1 . N 1 

' - A = _^ = expiy-(2Tr(p-r-Tr(Q-f) I (2.21) 



det(i-P-) [^1^ 



Notice that the kernels P and Q have different support. In can be checked order by order, 
that the above formula results in 



i)ip{92-e3)...ip{en-ei) (2.22) 



Adding this constant to (2.18) one finds 

\og{ga) = —j d9Qaa{e)\og{l + e~'^''^) + -\ogA (2.23) 

where 

Qaa{9) = ^a{e) " ^{26) - 7r5{9), MO) = -i^ log Ra{e) 

This is the exact non-perturbative g'-function as obtained in [14]. We wish to note that the 
equivalence of (2.21) and (2.22) was already observed in [15]. 

This derivation does not take into account the fluctuations around the saddle point, which 
are expected to give a non-extensive contribution [17]. However, equation (2.16) already looks 
like a periodic be. quantization (with an additional factor) and it is generally accepted, that there 
is no 0(1) piece in the periodic case. This implies that eq. (2.23) already takes into account all 
0(1) contributions, ie. there are no additional terms due to the fluctuations. 

The reader may be concerned that our arguments are not convincing and that the proposed 
re-normalization might be dictated by the result itself. We agree that our derivation of (2.23) is 
certainly not rigorous, however, it already shows that the density of states play an important role 
in the definition of the partition function. In the next section we present a general calculation 
built on first principles, which properly takes into account both the density of states and the 
fluctuations around the saddle point. In the massive case we obtain the same result (2.23). 



3 0{l) pieces to the free energy — formal derivation 

In this section we present a general framework to determine the 0(1) pieces to the free energy. 
We consider Bethe-Ansatz systems with one particle type and one rapidity variable 9 which takes 
its values from a domain B. The Bethe-Yang equations are assumed to take the form 

^iQ,ie„...M ^ -Q 5By(e„ 0^) = 1 (3.1) 

where pj = p{9j). The pure number a and the phases RBY{(^j) and SBY{dj,0k) depend on 
the problem at hand. In massive relativistic models (or arbitrary non-relativistic models) with 
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periodic boundary conditions we have 



a = l Rby{0) = 1 

SBY{(^ji(^k) = S{6j - 9k) 
In the case of integrable integrable boundaries 

a = 2 RBYi9) = Ra{0)Rb{0) 
SBY{9j,9k) = S{9j — 9k)S{9j + 9k) 

For the treatment of massless relativistic models see section 5. It is assumed that the unitarity 
condition is satisfied in the form 

SBY{9j,9k)SBY{—9j,—9k) = 1 

To investigate the thermodynamic hmit, we spht up the 9 axis into intervals A9 and we use the 
densities Po{9) and ph{9) for the occupied states and the holes. The total density of states is 
p{9) = Po{9) + Ph{(^)- It follows from (3.1) that the densities satisfy the constraint 

1 f rl9' 

Ph{e) + Po{9) = ^^{0) + J^—K,{9,9')p^"H9'), (3.2) 

where 

Ki{9,9') = -i^logSBY{9,9') 

and 

a{9) = a-^p{9) + e{9) (3.3) 

with 

e{9) = -i-^^log(^RBY{9)SBY{-9,-9)) (3.4) 

Depending on the problem at hand there can be an extra —2tt5{9) term to (3.4) to take into 
account the Pauli-principle. 

Based on the calculations of the previous section we assume that the density of states for an 
A^-particle state can be written as 

N 

Pn{9i, ...,9n)=MxII {2TrLp{9,)) (3.5) 
i=i 

where is a bounded finite number for a smooth distribution of roots. In the thermodynamic 
limit it is given typically by a Fredholm-determinant; the explicit calculation for a generic case 
will be presented later in this section. 

The partition function Z is obtained as a summation over all possible momentum quantum 
numbers in a finite volume L. In the thermodynamic limit it is evaluated by a functional integral 
over all possible density functions. The usual prescription to define a regularized functional 
integral for Z is 

E ^ I ...lll{LA9dpo{9)) (3.6) 

configurations 6 

Here the product is over the discretization points and there is one integral for the density at that 
point. 
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We believe that this prescription is not correct for Bethe Ansatz systems. The partition 
function depends on the configuration which minimizes the free energy and also on the number 
of available states for that particular configuration. The calculations which result in formulas 
like (3.5) show that the number of available states for a Bethe Ansatz system does not behave 
like the density in a free theory. Instead, there is the non-trivial finite factor M which has to be 
taken into account. In other words, the constraints for the configuration space of the theory have 
to be incorporated in the definition of the functional integral. 

In the usual evaluation schemes the entropy considerations are worked out in rapidity space, 
and the functional integral is performed using (3.6). In this procedure one is free to choose the 
discretization mesh A6 in a fairly wide range (for a discussion on this point see [4, 17]). The 
calculations are safe in the sense that the end result does not depend on the actual value of A9. 
However, A9 should be chosen large enough so that there are enough particles in the interval 
{9, 9 + A9) to make the entropy considerations meaningful. On the other hand, the "ultimate 
discretization" would be to have one discretization point for each particle in the system; in this 
case the variation of the free energy would be performed on the level of the momentum quan- 
tum numbers. Such a scheme is obviously not amenable for practical purposes. The solution is 
to perform a transformation from the space of momentum quantum numbers to the space of 
rapidities; the important point is to keep the Jacobian associated to this mapping. In the ther- 
modynamic limit the Jacobian behaves as given by (3.5), where the product over j = 1 . . . N can 
be interpreted as the "ultimate discretization" of the functional integral (3.6). A very important 
observation is that in the thermodynamic limit the pre-factor Af does not depend on the number 
of particles involved, therefore we conclude that it must be present in any discretization scheme. 

Based on the above considerations we propose the following definition of the functional inte- 
gral: 

j ... j MX[[LA9dpo{9)) (3.7) 

configurations 9 

The factor N depends on the particular thermodynamic configuration which minimizes the free 
energy functional. However, experience shows that 7\A is a finite number which does not grow 
with the volume. Therefore it is expected that it will not shift the position of the saddle point 
and its only effect is to correctly normalize the partition function. It follows that the prescription 
(3.7) is equivalent to 

M j ... j\[{LA9dpo{9)) (3.8) 

configurations 8 

and the difference between (3.7) and (3.8) vanishes in the thermodynamic limit. 

In the following we evaluate the partition function using the prescription (3.8). The calculation 
proceeds in four steps: 

1. Establishing the Thermodynamic Bethe Ansatz using the usual prescription (3.6) 

2. Evaluating the contribution of the quadratic fluctuations around the saddle point 

3. Evaluating the normalization factor M 

4. Writing down the final result with the correct normalization (3.8) 

The treatment of the quadratic fluctuations follows exactly the same way as in [17]; we simply 
restate the results of [17] in our general framework. In order the keep the exposition simple we 
omit some of the technical details and the checks of the various approximations we make. For a 
careful treatment the reader is referred to the original work ]17]. 
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step 1 

We evaluate the partition function using the functional integral (3.6); the entropy consider- 
ations are worked out in rapidity space. The number of micro- canonical configurations is given 

by 

e 

where 

fLA9p{9)\_ {LAep{e))\ 



UJ 



\LAepo{e)J {LAOpoie)) i {LAOphie)) \ 

We approximate it using Stirling-formula as 

\oguj = LAes{p{e)) + c;{p{e)) + . . . , 

where 

s{p{e)) = p{e) \og{p{9)) - pM \og{po{e)) - pm \og{ph{e)) 

and 

?(pW) = -^log(2vrLA0) + ^log- ^^^^ 



2 ' ' 2 ^ po{d)ph{e) 

The free energy functional can be written as 

F[p{e)] = lY, {e{e)po{e) - Ts{e)) ao (3.9) 

e 

The usual minimalization procedure yields the integral equation 

e{e)/T = e{e)+ f ^ Ki{e,9')log{l + e~<''^), (3.10) 
where we introduced the pseudo-energy function as 

e-'^'^ (3.11) 



Po{()) -e{9) 



Ph{e) 

The "minimal part" of the free energy can be expressed simply as 

F^^ = -LT I ^ a(0)log(l + e-^W), (3.12) 
Jb 

where cr{9) is given by (3.3). 
Step 2 

We consider the fluctuations around the saddle point solution. As a first step we expand the 
free energy functional (3.9) around Po{k) and ph{k): 



2\ po{e) + ph{9) ph{9) Po{9) I 

Here the quantities 



ro{9) = 6po{9) and rh{9) = 6ph{9) (3.14) 
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are constrained by 

ro{e) + rh{9) = J2 Ki {O, 9') ro{9')M' . (3.15) 
The absolute magnitude of the partition function is then evaluated as 

Z = M X e"^-^-- X [ ■■■ fYl{LA9dr{9))x (3.16) 



^\ ^[ 2 [ po{9) + ph{9) ph{9) Po{9)J 2 poPhLA927r J 
We introduce new integration variables ^(0) according to 



This leads to 

where 5g qi is the Kronecker symbol, and 



_ 1 po{9)p,,{9)Ae po{9')p,{e')Ae' 

- 7J9)i Po{9) + p,{9)'''^'^ ' Po{e') + pn{e') ^'-''^ 

Finally, changing the integration variable to ^ in (3.16) we have 



(3.20) 



= AAe-^^-'"(det[V-^e,e']) ' 
with 



Taking the limit A9 —?■ d9 yields 



-1 



Z = ATe"^-^™ ( det (i - Ki) ) (3.22) 



Step 3 

The density of states is given by the determinant 

dQi 



Pn{9i,- ■■ ,9isi) = det J, Ji 



ik 



09, 
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A straightforward calculation shows that the matrix elements are given by 

N 

i=i 

where 

K2{e,e') = i—\ogSBY{e,e') 

The matrix can be written as 

J' = GQ, where now 



with 

TV 



jj=La{6j) + J2Ki{ei,02) (3.23) 

i=l 

It follows from (3.2) that in the thermodynamic limit 

7j — ^ 2TrLp{6j) 
The elements of Gat can be written asymptotically as 

1 K2{9i,ej) 



Using (3.11) we conclude that the limit of detGjj is given by the Fredholm determinant with 
kernel K2{0, 6'). Therefore the normalization constant of the free energy functional (3.8) is given 

by 

= det (i - i^a) (3.24) 

Step 4 — Our main result 

Substituting (3.24) into (3.22) yields 

r.p det (i - Ks) , , 

^^g-/3Fmin L ^ (3.25 

det (l -Ki) 

Here Ki and K2 are integral operators which act on functions defined on B as 
and the kernels are given by 

Kiie,e') = -i^iogSBY{o,e') K2{e,e') = i-^\ogSBY{e,e') 

The phase shift Sby{0,0') is defined implicitly by the Bethe-Yang equations (3.1) and -Fmin is 
given by formula (3.12). 
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4 Explicit examples massive relativistic models 

In this section we evaluate (3.25) explicitly in four different cases; we restrict ourselves 
relativistic models. Non-relativistic ones can be treated in the same manner. 

4.1 Free fermionic gas 

In this case 

,B = M and a{9) = mcosh.6. 
There is no interaction between the particles, therefore 

Ki{e,e') = K2{e,e') = and M=i. 

The thermodynamics is trivial: 

£{e) = e{e)/T = mcoshe/T 

and 

^mcosh01og(l + e-"^^°^^''/^) 

-oo 

There is no 0(1) piece. 

4.2 Interacting particles with periodic boundary conditions 

In this case 

B = R and a {9) = in cosh 9. 

The scattering phase shift is 

Sby{9,9') = S{9-9'). 
Therefore the two integral kernels are given by 

Ki{9,9') = K2{9,9') = v{9-9') 

The free energy reads 

/oo Jfl 
^ ?n cosh log (l + e-^(^)) (4 
-oo 

Equation (3.25) results in 

det{i-K2) ^.pF". -^F". 

/j — g " min ^ g mm 

det (i - kx) 

The two Fredholm-determinants coincide and there is no 0(1) piece. 
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4.3 Periodic be. with a purely transmitting defect 

One can consider the Betlie Ansatz equations in the presence of a purely transmitting defect: 

where T(0) describes the scattering between the particles and the defect. Similar to the previous 
case one has 

However, the function (t(0) is now given by 

(t(0) = m cosh ^ + 99^(0), where = — i — — log T(0) 

L du 

The partition function is expressed as 

^ ^ det (1 - K2) ^-pfI^^-pft ^ ^-PFl,^-liFT 
det (i - ki) 

where Ft is an 0{1) piece given by 

/OO 
— (^TWlog(l + e-^(^)) 
-00 ^'^ 

This is in agreement with the previous results in the literature (see for example [23]). 

4.4 Integrable boundaries 

Here the domain of integrations is i3 = M'^ and the scattering phase shift is given by 

SBY{o,e') = s{6 -9')s{9 + e') 

The integral kernels are given by 

Ki{6,6') = 'f{9-6') + ^{e + 6') 

K2ie,e') = ip{e-e')-ip{e + e') 

The associated Fredholm-determinants are det (l — Q^) where the operators are defined as 

(Q±(/))(x) = - y) ± ^(^ + y)) (4-2) 

The function (t{9) is given by 

a{e) = 2mcos\ie + Qab{0) 

where 

Qam = -i-^^logRab{e) - 2^5{e) 
The TEA equation (3.10) takes the form 

(•CO jal 



m 



roo jal 

cosh e/T = e{0) + / ^ {^{9 - 0') + ip{9 + 6')) log (l + e'^^^')) 
Jo 27r 
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One can define the pseudo-energy also for negative values of rapidities as e{9) = e{—9), then the 
above equation above is equivalent to 

/oo jg' 
-oo 2vr 

which is the usual periodic boundary conditions TBA. The minimal part of the free energy is 
expressed as 

^■OO Jfl 1 

F^in = -LT — (2m cosh 9 + -6,6(0)) log (l + e'^^^)) = F^,^ + Fab 
where F^^^ is the free energy of the periodic be system (4.1) and Fab is an 0{1) piece given by 

/OO Jfl 
— e,fe(0)log(l + e-^W) (4.3) 

The partition function reads 

det (i - Q+) 

Therefore the 0(1) piece to the free energy is given by 

log(5a5fe) = -liFab + log '1''^ " (4.4) 

det (1 — Q+j 

In Appendix A we show that the above ratio of Fredholm-determinants reproduces the constant 
log^ defined in (2.22). Choosing the boundary conditions a = b the final result can be written 

as 

1 /"OO JQ 

log(9.) = 2 y_ log {1 + e-'<'y) + (4.5) 

+5 E i £ t . . . £ t (n rr^) ^(«. + ^^M^^ - «3) ■ ■ ■ - «o 

This is the exact non-perturbative (^-function as obtained in [14]. 



5 Explicit examples — massless relativistic models 

In this section we consider massless relativistic theories with diagonal scattering [24, 25]; 
for an introduction to massless scattering the reader is referred to [26]. We assume that there 
is only one particle type in the spectrum. There is an effective doubling because one has to 
treat left-moving and right-moving particles separately. We consider parity-invariant theories; 
the S-matrices are given by 

Sll{9) = Srr{9) Slr{9) = Srl{9) 

The TBA equations are usually written down for periodic boundary conditions, in which case 
the TBA is a two-component system: 

em = ^rnRe' -J ^ [^u{9 - 9') log (l + e'^^^^')) + ^12(0 - 9') log (l + e'^^^^'))) 

1 /"OO Jfjl , ^ 

e2{9) = -mRe~' -J — (^^2i{9 - 9') log (l + e^^^^^')) + (^22(0 - ^0 log (l + e^^^^'''))) 

(5.1) 



19 



The kernels are defined as 



Vim = ip22{e) = -i^iogSRRie) ipi2{9) = V2i{e) = -i-^iogSLde) 

One can then use the symmetry ei(^) = E2{—0) to transform the above system into a single 
equation: 

e{e) = -mRe'- / ^{^^^{9 - 9') + vi2{e + e'))\og{l + e~'^''^) (5.2) 

^ J~oo ^'^ 

The free energy is expressed as 

^|^me^log(l + e-W) (5.3) 

It is generally expected that for periodic boundary conditions there is no 0(1) piece. 

In the following we apply the formalism of section 3 to determine the (^-function in the 
presence of two integrable boundaries with reflection factors Ra{0) and Rb{0). Contrary to the 
periodic case, in the open system there is no distinction between the left-movers and the right- 
movers. When a left-mover scatters off the left-boundary, it becomes a right-mover with the same 
energy and reversed momentum. Therefore both the Bethe Ansatz and also the thermodynamics 
can be written down in terms of only one particle species and one pseudo-energy function. 

We use simple heuristic arguments to write down the Bethe- Yang equations. The particles 
will be parametrized with the rapidity variable ^ G M which refers to the situation when the 
particle is moving to the right, ie. 

1 1 

e{e) = -me^ p{e) = -?ne^ 

When a particle with rapidity 6 is taken back and forth in the open system it meets every other 
particle twice. The two scattering processes are described by Srr{9 — 6') and Slr{0 + 9'). The 
Bethe- Yang equations read 

e'^P^'^Ra{e)Rb{e) n SRR{9k - 9j)SLR{ek + %) = 1 (5.4) 

To establish the connection with our general formalism we extract the quantities 

a = 2, Rby{0) = Ra{e)Rb{e)SLR{-2e), 

and for the scattering phase shift we find 

Sby{0, 9') = SRR{e - 9')Slr{9 + 9') 

which yields the integration kernels 

Ki{9,9') = vu{9 - 9') + vi2{9 + 9') 
K2i9,9') = vn{9 - 9') - vi2{9 + 9') 

The function (t{9) is given by 

a{9) = me^ + Qabi9) 
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where 

QabiO) = -i^J log Rab{9), Rab{0) = Ra{e)Rb{e)SLR{-2e) 

Note that contrary to (2.5) the 6{6) term is missing from cr{9) because there are no formal 
solutions with 6 = which should be canceled. 
Formula (3.10) yields the TBA equation 

e{9) = -mRe'- / ^{cpu{9 - 9') + ipi2{e + 6')) log {l + e-<''^) (5.5) 

This equation coincides with (5.2), although it was derived from a conceptually different Bethe 
Ansatz. The interpretation is straightforward: the distribution of roots (and therefore the 0{L) 
pieces of the free energy) do not depend on the boundary conditions, as it is expected on general 
grounds. 

Formula (3.12) yields 

Here -F^j^ is given by (5.3) and Fab is an 0(1) piece given by 

|^e,6(0)log(l + e"^W) (5.6) 

Notice the factor of 2 as compared to (4.3). 

Putting everything together, equation (3.25) yields 

det (l — K2) opp op , , 

Z = ) ^e-^^-i""''^'"' (5.7) 

det(l-Ki) 

In the next two subsections we explicitly work out the details for simple scattering theories 
with one particle species. The generalization to other models with more than one particles (for 
example the scattering theory in [27]) can be treated with the straightforward extension of (5.7). 
Finally we mention that the boundary independent part of (5.7) can be written in the form 



1 det (1 - K2) ^1/ „ 

log- — J- TTT = >^-(Tr(i^i - Tr A2, 

det(l-Ki) ^,nV 



The second summation runs over = 1,2, the pseudo-energies are given by 

£^{e) = e{9) e2{9) = e{-9) 

with e{9) being the solution of (5.5) and the kernels are defined as 

+ fo\ ( ^pi2{9) ior j = k _ ( ipii{e) for j = k 

^^^^'^ = [ for , / k ^d'^ = [ ^M9) for j + k 

Equation (5.8) can be proven term by term using the symmetry Lp^j^{9) = Lp^^{—6). 
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5.1 The massless flow from tri-critical Ising to critical Ising 

The simplest non-trivial massless model is probably the scattering theory describing the flow 
from the tri-critical Ising to the critical Ising model [24] . In this theory there is only one particle 
species and the scattering is described by 

SLL{e) = SRR{e) = 1, SLR{e) = - tanh(0/2 - m/4) (5.9) 

In [16] ^-function flows were studied between different conformal boundary conditions of the UV 
and IR theories in those cases where both the bulk and the boundary perturbations are integrable 
and compatible with each other. There are two such possibilities: 

1. The flow from the boundary condition (0+) of tri-critical Ising to (-1-) in Ising 

2. The flow from the boundary condition [d) of tri-critical Ising to (/) (free) in Ising 

Both flows are induced by the $13 perturbation on the boundary. For the precise definition of 
the boundary conditions we refer to [16] and references therein. Here we show that the results of 
[16] can be derived from our general formalism. Most importantly, we present an all-orders proof 
of the boundary-independent part of the ^(-function, which differs from the massive version. 
In the present case the integration kernels are given by 

Ki{9,e') = -K2{e,9') = ^{6 + 9')^ 
a{e) = me^ + Qab{0) 



' cosh(e + Q') 
The function o"(^) reads 



where 

Qah{Q) = log RahiQ). Rab{0) = Ra{0)Rb{0)SLR{-2e) 

The possible reflection factors Raid) and Rb{9) were specified in [16]. 
The TBA equation is given by 

e{9) = -mRe^- —^i9 + 9')log{l + e-'^^'^) (5.10) 
2 J —CO 27r 

Formula (3.12) yields 

Here F^^^ is given by (5.3) and Fab is an 0(1) piece given by 

/OO JO 
£e,,(0)log(l + e-(^)) 

Putting everything together, equation (3.25) yields 

det(l-FP+) npp np 

Z = \ ^e-'^^---'^^"'' (5.11 

det (1 - P+) 

where the operator P+ acts on functions defined on M as 
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The ratio of Fredholm-determinants can be evaluated using (2.20): 



log ^!ili±i:^ = f 1 (TriP^f _ Tr(-P+' 
det(l-P+) ^1™^ 



We find the exact (^-function 

log g = log ga + log 50 
where the boundary dependent part is 



-(v,,(0)-<^(20))log(l + e--(^)) (5.13) 



with 

¥'a(0) = -i^logiia(^) 

and the boundary independent part is 



log 50= (5.14) 

I ih ^ ■ ■ ■ ^ ( n iriiso) - -(^-' + 

Equation (5.14) is in agreement with the corresponding formula of [16]. However, (5.13) coincides 
with the result of [16] only in the case of the flow to the (/) boundary condition in the Ising 
model. In the other case, namely the flow from (0+) in tri-critical Ising to (+) in Ising there is 
a missing term — ^ log 2. This discrepancy can be explained by the fact, that in the IR limit the 
boundary condition corresponds to a microscopic theory where the ground state degeneracy of 
2 is removed. The scattering theory describes the variation of the ^-function with respect to the 
temperature, therefore it is natural to assume that the extra term has to be added not just in 
the IR limit (which corresponds to zero temperature), but also for the whole RG flow. Therefore 
we write 

log5a = -2log2 + y_ — (v.„(0)-v.(20))log(l + e-^W) 

which is in agreement with [16]. 

Finally we mention that (5.14) follows from the general formula (5.8) after substituting 

¥'ii(e) = (^12(0) = V3(e) 

and making the appropriate change of variables. 



5.2 The massless flow Ms^s + $2,1 M2,5 

In [28] a simple massless scattering theory with one particle species was proposed to describe 
the flow from the minimal model M3 5 to Mi^^ induced by the perturbing field ^>2,i [29]. In this 
model the scattering is described by 

= SRR{e) = Sly{0), SLR{e) = SRUe) = {Sly{0))-\ (5.15) 
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where Si,y{0) is the S-matrix of the massive Lee- Yang model [30, 31] 



^ sinh ^ + i sin(7r/3) 

sinh 9 — i sin(7r/3) 

The massive Lee- Yang model is the <I>i^3 perturbation of the minimal model M2.5. In massless 
theories the LL and RR scattering matrices are scale-invariant and they describe the IR limiting 
CFT; this was the motivation for the choice of Sll and Shr in (5.15). 

The possible reflection factors of this model have not yet been written down. Nevertheless it 
is useful to derive the g-function, leaving the factors Ra{0) and Rb{0) unspecified. The bound- 
ary dependent part will be given by (5.6); in the following we concentrate on the boundary 
independent part. 

Given the scattering matrices (5.15) the integral kernels are given by 

Ki{6,e') = ^{d-e')-^{e + 6') 
K2{e,e') = ip{e-e') + ^{6 + e'), 



where 



The TBA equation reads 

e{9) = -mRe^ - — {^{9 - 9') - ^{9 + 9')) log (l + e'^^^')) (5.16) 
One can use the general formula (5.8) to express the boundary independent part of the (7-function 



as 



1 ^ det (1 - K2) 
log 50 = 7: log 



2 det (1 - A'l) 

r°° rift fJL 1 \ 

ip{9i + 92M92 - ^3) • • • HOn 



The summations run over aj = 1,2 and the pseudo-energies are given by 

£^{9)=e{9) £2{9) = e{-9) 

where e{9) is the solution of (5.16). Notice that due to the specific form of the S-matrix (5.15) 
there is an overall factor of (—1) as compared to the massive case (2.22). 



6 Conclusions 

We have studied the partition function of Bethe Ansatz solvable models as a function of the 
volume and the temperature. We have shown how to obtain 0{LP) pieces to the free energy 
in the framework of the Thermodynamic Bethe Ansatz: our main result is equation (3.25). In 
addition to possible boundary dependent parts incorporated in the "minimal part" of the free 
energy Fmin the formula (3.25) involves two Fredholm-determinants which depend only on the 
scattering in the bulk. In relativistic boundary field theory these two pieces are responsible for the 
boundary-independent part of the (^-function. We have presented a new result (5.7) which applies 



24 



to massless relativistic theories with arbitrary diagonal scattering in the bulk. This formula could 
be used to study massless bulk-boundary flows along the lines of [16]. 

Formula (3.25) can be applied in a very straightforward way once the Bethe-Yang equations 
have been established. Therefore it is very natural to conjecture that a similar result will hold in 
theories with non-diagonal scattering. In these models the finite volume quantization proceeds 
through the diagonalisation of the transfer matrix, which can be achieved by the introduction of 
the so-called magnonic (or spin) particles [2, 32, 33, 34]. Once this algebraic problem is solved, 
the derivation of the TBA follows straightforwardly; a common property is that there are no 
energy-terms e{9) associated to the magnonic modes. We believe that our arguments can be 
applied to these Bethe Ansatz systems, particularly in those cases when the TBA results in a 
finite set of equations. The study of these ^r-functions with magnonic modes is left for future 
work. 

Also, it would be interesting to study the flow of the excited states quantities G^{R) defined 
in (1.3). In models with discrete symmetries some of the excited states can be treated simply by 
introducing real or complex fugacities = e"^"/"^ for the different particle types [35, 36]. We 
expect that our present results will hold in these modified TBA systems. However, the problem 
of particle-like excited states will probably require new ideas. In [37] it was shown that generally 
the amplitudes {R) include normalization factors which are analytic in 1/R; it is not yet clear 
whether these factors can be interpreted in the framework of Thermodynamic Bethe Ansatz. 

Acknowledgements: We would like to thank L. Palla, Z. Bajnok, M. Kormos and in partic- 
ular R. Tateo and G. Takacs for encouraging and very helpful discussions. Also, we are grateful 
to M. Kormos, G. Takacs and G. Palacios for useful comments on the manuscript. 



Appendix 

A Relations between the Fredholm-determinants 

Here we consider simple relations between the different Fredholm-determinants introduced 
in the main text. 

The operators and P~ act on functions defined on M as 



The function ip{0) = ■d'{9) is the scattering kernel. A very important property is that it is 
symmetric: 



Let us decompose the real line as M = -|- M (the point x = has zero measure, therefore it 
is irrelevant in the present context). It is easy to see, that in this decomposition the operators 





(g± {g)) (x) = j^^ {^{^ - y) ± ^(^ + y)) j^^ff(y) 
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and P can be written in block form as 

where A and B are integral operators on with kernels {p{x — y) and ip{x + y), respectively. 
The determinant of the operators 1 — P^ can be evaluated as 

det (i - P+) = det (i - (i + B)) det (i + (i - S)) 

det (i - P~) = det (i - (i + i?)) det (i - (i - B)) 



Using the relations 



one gets 



Q+ + Q" ^ _ Q+ - 



det (1 - P+) = det (i - Q+) det (i + Q") (A.2) 
det (i - P-) = det (i - Q+) det (i - Q") (A.3) 

The relation (A.3) was used in the main text to prove the equivalence of the expressions (2.23) 
and (4.4). We wish to mention that eqs. (A.2)-(A.3) can be proven alternatively term by term 
using formula (2.20). 
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